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Abstract — We investigate power allocation for users in a spec- 
trum underlay cognitive network. Our objective is to find a power 
control scheme that allocates transmit power for both primary 
and secondary users so that the overall network throughput is 
maximized while maintaining the quality of service (QoS) of the 
primary users greater than a certain minimum limit. Since an 
optimum solution to our problem is computationally intractable, 
as the optimization problem is non-convex, we propose an 
iterative algorithm based on sequential geometric programming, 
that is proved to converge to at least a local optimum solution. 
We use the proposed algorithm to show how a spectrum underlay 
network would achieve higher throughput with secondary users 
operation than with primary users operating alone. Also, we show 
via simulations that the loss in primary throughput due to the 
admission of the secondary users is accompanied by a reduction 
in the total primary transmit power. Q 

I. Introduction 

Cognitive radios, or secondary users, are spectrum agile 
devices that can sense the spectrum, identify the underutilized 
bands, and opportunistically share those bands with the origi- 
nal users (primary users) such that they do not cause harmful 
interference on those primary users (TJ. Hence, more users can 
share the spectrum efficiently and provide more utilization to 
the underutilized bands. However, when a secondary user is 
severely interfering with the primary users, it should abandon 
this primary band and try to exploit another underutilized one. 

Cognitive radios can operate with a primary network using 
either overlay, or underlay spectrum sharing Q. In spectrum 
overlay, secondary users are allowed to access the primary 
bands that are totally free of primary operation. Consequently, 
the interference caused on primary users is almost negligible, 
and robust sensing is demanded of secondary users. In spec- 
trum underlay, secondary users can share the same spectrum 
band concurrently with the primary users, so long as the 
interference caused on the primary users is small enough to 
guarantee minimum primary QoS. The secondary users do not 
have to employ sophisticated sensing mechanisms. Instead, 
efficient resource allocation procedure is needed so that the 
spectrum can be shared with high utilization. 

In this paper we consider the problem of power allocation 
for primary and secondary users in a spectrum underlay sys- 
tem. Most of the existing work separates the power allocation 
for the secondary network from that of the primary network. 
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In [3 1 and [4] the authors consider the power allocation for 
secondary users under maximum tolerable interference con- 
straint on the primary users, besides minimum QoS constraint 
on each secondary user. They address two scenarios, the first 
is the maximization of the secondary network throughput in 
case that all of the secondary users can be supported under 
the given constraints. The second is the maximization of the 
number of admitted secondary users under their minimum 
QoS constraints when not all of them can be supported. 
The second scenario is also addressed in [5] where a lower 
complexity and more efficient algorithm is proposed. In O the 
authors combine the power allocation and admission control 
for secondary users in one problem aiming at maximizing the 
overall secondary throughput. 

In this work we consider an underutilized primary network 
and investigate the overall network throughput in the presence 
of cognitive radios. Our objective is to jointly allocate power 
to both primary and secondary users so that the overall 
network throughput is maximized, while protecting the QoS 
of the primary users by keeping their signal-to-interference- 
and-noise-ratio (SINR) above a minimum level. Since power 
allocation for maximum sum throughput in multiuser wireless 
networks is a non-convex problem, we introduce an iterative 
algorithm that is proved to converge to an optimal value (not 
necessarily the global) of the total sum throughput. We inves- 
tigate the efficiency of the solutions to which the proposed 
algorithm converges via numerical simulations, and show that 
the network can achieve better utilization by incorporating 
secondary users operation. Also, we show the throughput and 
power prices paid by primary users in order to sustain the 
secondary operation and achieve higher utilization. Our results 
show that although the primary users may lose throughput 
due to the secondary operation, their total transmit power is 
significantly reduced. 

The rest of this paper is organized as follows. In Section 
HIl we describe the underlay network model. In Section [Til] 
we state the problem formulation. We introduce the proposed 
algorithm and establish its convergence in Section [IV] Sim- 
ulation results are presented in Section [V] and the work is 
concluded in Section [VI] 

II. System Model 

We consider a spectrum sharing model where secondary 
users share the same frequency band with primary users. In 



this model we focus on the uplink transmission from primary 
users toward their base station (BS), whereas the secondary 
users are communicating with each other in an ad-hoc fashion 
as shown in Fig. Q] 




primary user does not fall below a minimum level. The QoS 
constraint on the primary users are 
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This is equivalent to guaranteeing a minimum rate for each pri- 
mary user. Moreover, each transmitter, secondary or primary, 
cannot increase its transmit power beyond a certain maximum 
level. The power constraints on secondary and primary users 
are 
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Using the definition of the throughput of user i in bps/Hz 
as log 2 (1 + SINRi), the optimization problem is written as 
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Fig. 1. System Model 
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The primary network is assumed to have a single BS with 
K primary users. We denote the transmit power of the j th 
primary user by p P , j = 1, 2, • • • , K, the channel gain from 
user j to the base station by and the SINR of primary 
user j by x P . Thus, 
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where there exist N secondary links within the coverage area 
of the primary BS. Each secondary link i comprises one 
transmitter and one corresponding receiver. The parameter pf 
is the transmit power of the transmitter of secondary link i, 
<?Q S j is the channel gain between the i th secondary transmitter 
and the primary BS, and No is the background noise power 
in the frequency band of operation. The channel gain between 
a secondary transmitter i and a secondary receiver j is g^\, 
whereas the channel gain between a primary transmitter i and 
a secondary receiver j is g^ p . Hence, we can write the SINR 
of the secondary link i, x, 
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A central controller is assumed to be responsible for the 
resource allocation for the network nodes. This controller is 
assumed to have knowledge of all of the channels in the 
system. 

III. Problem Formulation 

For the given system, we consider the sum throughput. 
The network users are allocated powers so that the total 
throughput, primary and secondary, is maximized. However, 
in order to guarantee minimum QoS for the primary users, the 
power control scheme should guarantee that the SINR of each 
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where the optimization variables are p p and pf for j = 
l,--- ,K and i = , AT. The SINRs: x P and xf are 

defined in (fl]i and (fJJ. Parameters a s and a p are control 
binary variables to switch between three different scenarios: 
1) a s = 0, a p = 1 in which secondary users are not allowed 
to operate, 2) a s = 1, a p = 1 in which secondary users can 
share the spectrum with the primary users, and 3) a s = 1, 
a p = in which the secondary sum throughput is maximized 
under the primary QoS constraints. We consider the second 
scenario in the next section. The optimization problem in (O 
is a non-convex one and has no tractable solution Q. In 
|E1, a similar optimization problem is reduced to a convex 
problem, by omitting the unity term added to the SINR in each 
throughput term. That is, the system is assumed to operate in 
high SINR regime where log 2 (1 + SINR,) w log 2 (SINRi) 
and the problem is then converted to a geometric program 
and solved for a global optimum. But the difference in our 
formulation is that there are no constraints on the SINR of the 
secondary network, i.e., the system is opportunistic for the 
secondary users, consequently high SINR is not guaranteed 
for every user. 

Also, the feasibility of the problem in (O necessitates that 
the primary network alone is feasible. That is, there exists a 
power allocation for every primary transmitter so that all the 
primary users can be supported exactly at their minimum QoS 
requirements. We assume that the proposed algorithm takes 
place when the primary network is feasible. 

IV. Iterative Geometric Programming 

We convert the non-convex problem in © to a convex one 
by approximating the objective function by a monomial. We 



formulate the resulting problem as a geometric program that 
can be solved efficiently with polynomial computational com- 
plexity. Then, we iteratively update the approximate function 
and solve the geometric program again till convergence to at 
least a local maximum solution of the original problem. 

A. Geometric Programming 

Before defining the geometric program we first introduce 
the following two terms: 

> A monomial m(x) = Cx^x^ 2 ■ ■ ■ x® n where C is non- 
negative real constant, and a\,a,2,-" : a n are rea l con- 
stants. 



• A posynomial s(x) = X^fc-i CkX°^ k x. 
sum of monomials (8). 

A geometric program is defined as follows: 
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consequently, the approximate optimization problem will be 
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minimize /o(x) 

subject to /j(x) < 1, i = 1, 2, • ■ • , I (7) 
hj(x) = 1, j = 1, 2, • • • , m 

where, x S /i( x X i = 0, 1, • • * ,2 are posynomials and 

hj(x), j = 1, 2, • • ■ , m are monomials. 

Geometric programs can be simply converted into convex 

problems 0. Since Y,i=i lo S2 /fat) = lo S2 Il"=i /fat). the 
optimization problem (O can be written as: 
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where p s and p p are the vectors of secondary and primary 
power allocations, respectively, and x s and x p are vectors 
of SINR allocations of secondary and primary users. We 
note that the constraints of the optimization problem ([H) 
satisfy the conditions of posynomials. However, the conversion 
of the objective from maximization to minimization while 
maintaining the posynomial form is not achievable. For this 
reason we approximate the objective function, /o(x s ,x p ), to a 
monomial /o(x s ,X p ) where 
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where the parameters c, Af, and X P , i — 1, • • • ,N and 
j = 1 , • • • , K are determined iteratively so that, for any 
two successive iterations k and k — 1 we always have non- 
decreasing sum throughput, i.e., 

f (xf~ 1 \4 k -V) < f (xf\x p W) (11) 



where (Xq , ) is the solution to the approximate problem 
([Tol l for the k th iteration. 

B. Determination of fo(x s ,jc p ) for Non-decreasing Through- 
put 

In this subsection, we consider the complete determination 
of /o(x s ,x p ) in a certain iteration so that the inequality of 
( fTTT i holds. Let (x§,Xq) be the solution to the approximate 
problem ( TTOb in the previous iteration. Then, to have a sum 
throughput in the following iteration that is at least no less 
than that corresponding to (xq,Xq), the approximate function 
/o(x s ,x p ) should satisfy the following conditions: 

1) f Q (x s ,xP) > fo(x s ,xP) for all x s ,x p , 

2) / (xg,xg) = / (xg,xg), 

3) V/o(xg,xg) = V/ (xg,xg). 

The above three conditions are sufficient to guarantee that 
the solution of each approximate problem increases the ob- 
jective function, and after convergence of the series of the 
approximate problems, the Karush-Kuhn-Tucker conditions of 
the original problem in (JHJ will be satisfied [9], [10|. Let 
x = [xx,x 2 --- ,x N+K ] where Xj = x P , j = 1,2, ••■ ,K, 
xj = xf, j = K+i, i = 1, 2, • • • ,N, A = [Ai, A 2 , • • • , X N+K ] 
where Xj = X P , j = 1, 2, • • • , K, Xj — A|, j = K + i and 
' 1.2.---.A. 

From the third condition we have: 
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From the second condition we have: 
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By obtaining c, Af and for i = 1, • • • , N and j = 
1, • • • ,K, the approximate objective function /(xq,Xq) is 
determined for the next iteration. What remains is to show how 
would the new approximate function satisfy the first condition. 
It is required to prove that 

N+K N+K 
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Theorem 1: Fix i, for all > the following inequality is 
satisfied: 

l + ii>-75T-ii j Xi > (15) 

Proof: By combining both sides of ( fTBI ). and taking the 
log we define: 
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Now we prove that g(xi) < by showing that the maximum 
of g(xi) is zero. First, the stationary point is obtained by: 
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By substituting with x* in ( fTol i: 

g(x* = x Qi ) = 

The second derivative of g(xi) is: 

d 2 g(x*) _ -l/xoi 
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which is negative for XQi > 0. Therefore, g(xi) < for Xi > 
0, and consequently the inequality of $15[ is satisfied. 

Since the inequality of (TT3T > is satisfied for all i — 
1, 2, • • • ,N + K, inequality (fT4l > is also satisfied, and the first 
condition on the approximate function is maintained. ■ 

Our proposed iterative algorithm is implemented as follows: 

1) Choose positive initial values for c,Ai, i — 1, • • • ,N + 
K. 

2) Solve the approximate problem in (tlOt . 



3) If a target accuracy is reached, stop the algorithm. 
Otherwise go to step 4). 

4) Calculate new values of and c from ( fT2l . ( fT3l . then 
go to step 2). 

Note that, the proposed technique can be applied to any 
sum throughput maximization problem where the constraints 
are convex. In ifTTl . an approach to solve the sum throughput 
maximization in interference channels is proposed based on 
single and double condensation methods (51 • The single con- 
densation method is shown to be outperformed by the dou- 
ble condensation method. However the approximate objective 
function of the double condensation method was not proved 
to satisfy the first condition in Section. [IV] 

V. Numerical Results 

The simulation setup we use is close to that used in |4|. We 
assume that the secondary transmitters and the primary users 
are located in an area of size 1000m x 1000m with BS of the 
primary network located at the center. The receiving node of 
each secondary link is placed randomly in a 500m x 500m 
square with its transmitting node at the center. The channel 
gain between any transmitter j and receiver i is modeled as 
gij = Kq x lO^'j/ 10 x d^j, where dij is the corresponding 
distance, fiij is a random Gaussian variable with zero mean 
and a standard deviation of 6 dB to account for shadowing 
effects, and Ko — 10 3 is a factor that includes some system 
parameters such as antenna gain and carrier frequency. The 
background noise power is iVo = 10~ 12 watt. The maximum 
transmit power on each secondary transmitter is p smax = 0.1 
watt, and on each primary transmitter pP max — 1 watt. For 
simulations we use the CVX, a package for specifying and 
solving convex programs [12]. 

A. Convergence of the Iterative Algorithm 

We compare the value of the sum throughput to which 
the proposed algorithm converges with the optimum solution 
obtained by exhaustive search, for this simulation we consider 
N ~ 2 and K = 2, the minimum QoS for the primary 
transmitters is chosen to be j p = 6 dB. For the feasibility of 
the primary network we assume primary and secondary users 
employ CDMA with primary processing gain PG P = 20, 
and secondary processing gain PG S = 80. We run five 
different simulations where in each one the locations of all 
users are generated randomly. The exhaustive search result 
is obtained by discretizing the transmit power range of each 
user into 100 samples and trying all possible combinations 
of these samples, then picking the one that achieves the 
maximum sum throughput under the QoS constraints. We start 
the iterative algorithm with all exponents in / equal to one, 
i.e., A = [1, 1, 1, 1]. Fig. [2] shows how the iterative algorithm 
almost converges to the optimum solution for small number 
of users. The proposed algorithm, however, is not guaranteed 
to converge to the global solution in general situations. The 
convergence in a low SINR regime is shown in Fig. [3] where a 
simulation is run under unity processing gain for both primary 
and secondary users. And the minimum QoS for the primary 
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Fig. 2. Total network throughput in different simulations for both the iterative 
and exhaustive search schemes N = 2 and K = 2 
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Fig. 4. Total underlay network throughput for N = 3 and K - 
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Fig. 3. Convergence of the proposed algorithm in low SINR regime for 
N = 2 and K = 2 
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Fig. 5. Primary network sum throughput for N = 3 and K = 5 



users is j p = — 3 dB, while the initial condition is not changed. 
Here, we note that, the initial condition of equal values of A 
means that the system is at high SINR which is not true in the 
current scenario, and that is obvious from the solution of the 
first iteration, however the proposed algorithm then iteratively 
increases the throughput till it reaches a local maximum that 
is the exhaustive search solution in this setup. 

B. Network Utilization and Primary Cost 

In this subsection we compare between the total throughput 
of the primary network alone, and the total throughput of 
the overall underlay network. In this simulation we assume 
N = 3 users, and K = 5. Both primary and secondary 
users use processing gain of 80. The minimum QoS for each 
primary user is ^ p = 8 dB. We use the proposed iterative 
algorithm in the three scenarios discussed in Section [Ell] The 
initial components of the vector A are equal. We run five 
experiments, where in each one, the locations of all users 
are generated randomly and the three optimization scenarios 
take place. For each optimization scenario we calculate the 
total sum throughput, the primary sum throughput and the 



total primary transmit power. The results showing those three 
quantities are depicted in Fig. [4] Fig.|5]and Fig. [6] respectively. 
From Fig. [4] it is clear that the second scenario in which the 
primary-plus-secondary sum throughput is maximized is the 
one that leads to higher utilization in terms of total underlay 
network throughput. In Fig.|5]the primary throughput is shown, 
and it is clear that the concurrent operation of the secondary 
and primary users leads to degradation of the throughput 
that the primary network can achieve with no secondary 
operation. However, the throughput loss of the primary users 
in the second scenario is not severe and a cooperative primary 
network may tolerate this throughput price for the sake of 
more utilization. Fig. [6] shows the total transmit power of 
the primary network in each of the three scenarios. Although 
the primary users lose some throughput due to the secondary 
operation, their total transmit power is reduced significantly as 
a consequence of the throughput loss. Since the third scenario 
yields minimum transmit primary power under maximization 
of the secondary network throughput, the third scenario can be 
used in an optimization problem that aims at maximizing the 
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Fig. 6. Total primary transmit power for 7Y = 3 and K = 5 



secondary sum throughput while the primary users are kept at 
their minimum possible transmit power with QoS constraints. 

VI. Conclusion 

We have considered the resource allocation for spectrum 
underlay cognitive radio networks. Our objective is to allo- 
cate power jointly to primary and secondary users so that 
the total underlay network throughput is maximized. In the 
optimization problem the primary users are protected against 
the opportunistic behavior by setting a minimum QoS con- 
straint for each of them. Since the objective function of the 
optimization problem is non-convex, we introduce an iterative 
algorithm based on successive solutions to approximate ge- 
ometric programs. We provide a mathematical proof of the 
convergence of the proposed algorithm and show that the 
algorithm converges always to one of the local optima of the 
original problem. Then we show via simulations the efficient 
results to which the algorithm converges compared to the 
exhaustive search method. Finally, we show how the algorithm 
is used to increase the overall network throughput (utilization) 
with relatively small loss in primary throughput and with 
considerable decrease in total primary transmit power. 
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